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Abstract 

We study the classical Antonov problem (of retrieving the statistical equilibrium 
properties of a self-gravitating gas of classical particles obeying Boltzmann statistics 
in space and confined in a spherical box) for a rotating system. It is shown that 
a critical angular momentum Ac (or, in the canonical language, a critical angular 
velocity tOc) exists, such that for A < Ac the system's behaviour is qualitatively 
similar to that of a non-rotating gas, with a high energy disordered phase and a 
low energy collapsed phase ending with Antonov's limit, below which there is no 
equilibrium state. For A > Ac, instead, the low-energy phase is characterized by the 
formation of two dense clusters (a "binary star"). Remarkably, no Antonov limit 
is found for A > Ac. The thermodynamics of the system (phase diagram, caloric 
curves, local stability) is analyzed and compared with the recently-obtained picture 
emerging from a different type of statistics which forbids particle overlapping. 



1 Introduction 

The name "Antonov problem" is usually referred to the study of the mean- 
field equilibrium state of self-gravitating A^-body systems of classical particles 
in a finite volume. One considers the system with Hamiltonian 

1 ^ 

^ H^{{r,}, {p,}) = -Yp" + H{r^}) (1) 
^ 1=1 

l<i<j<Af I' * ' 3\ 

enclosed in a finite three-dimensional spherical box of radius R, which is nec- 
essary in order to prevent particles from evaporating. In (1,2): rj G F C 
and Pj e denote, respectively, the position and momentum of the i-th par- 
ticle, G is the gravitational constant, and V stands for the (volume of the) 
box containing the particles; we set the particles' masses to one. The main 
problems are: (a) find the equilibrium configurations, i.e. the particle density 
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profiles p(r) {fy p{r)dr — N) that maximize the microcanonical entropy {a — 
constant) 



{Dx — nili dxi)-^ and (b) analyze the corresponding thermodynamics. Clearly, 
the existence of such profiles and their actual structure will depend on the 
value of the total energy E. Intuitively, when E is large and positive (that 
is when the kinetic term dominates) a gaseous-type of structure will be pre- 
ferred, with particles uniformly spread over V , whereas when E is sufficiently 
negative particles are expected to attract each other in a dense core embed- 
ded in a vapor. The situation presents however several subtleties. Assuming 
that particles are point masses obeying Boltzmann statistics in space, i.e. that 
they can come arbitrarily close to each other, it is found within the mean-field 
approximation [1-5] that 

(a) there is no global entropy maximum, therefore all maxima are local and 
correspond to metastable situations; 

(b) if = e > ec — —0.335 (i.e. at sufficiently high energies, for we take R 
to be fixed) there exist local entropy extrema where the particle density 
satisfies the condition 



where A is a constant. Notice that (3 is itself a function of p. The physically 
admissible solutions of (4) are spherically symmetric (in fact, they can be 
mapped onto the solutions of the Emden equation y" + 2y' /x = e~^, with 
y(0) = y'{0) = 0) and are called microcanonical "isothermal spheres"; 

(c) in general, local extrema in the space of p's can be maxima (metastable 
states) or saddle points (unstable). Microcanonical isothermal spheres are 
local maxima if the "density contrast" p(0) / p{R) ^ 709, otherwise they 
are at most saddle points of the entropy surface; 

(d) if e lies in the range between —0.2 and Cc isothermal spheres have nega- 
tive specific heat ^ : in such states the system heats up and shrinks upon 
decreasing energy; 

(e) finally, if e < ec (i.e. at sufficiently low energies) there are not even local 
entropy extrema. 

The standard interpretation [2] is that, while at high energies an equilibrium 
configuration is achievable in the form of an isothermal sphere, below e,. the 
system collapses to a dense core and overheats, following the negative-specific- 
heat trend to its extreme consequences. Such a "transition" , occurring at the 

^ As first pointed out by Thirring [6], in the microcanonical ensemble self- 
gravitating systems can have negative specific heat, at odds with the canonical 
setup, in which it is related to energy fluctuations and positive definite. 




(3) 





3/2 



(4) 
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Antonov point e = ec, is known as "gravothermal collapse" or "catastrophe" 
(see also [7, 8] for a more formal treatment of this transition) . 

In this work we extend the above theory to rotating systems, by including 
the effects of the angular momentum L = Y^iLi T"? X Pj as a second conserved 
quantity besides E, and discuss the thermodynamics and the corresponding 
equilibrium states (which will now depend on both E and L). Our work is 
motivated mainly by the crucial role rotation is expected to play in astro- 
physical systems on different scales such as globular clusters and the cores of 
elliptic galaxies [9,10]. both of which can be modeled by (1,2). (Hence our 
"particles" can be interpreted equally well as atoms and as stars.) A gravity- 
driven collapse can be seen as the prime cause for the equilibrium shaping of 
self-gravitating systems, both rotating and not. But the nature and the result 
of this collapse are expected to change significantly in the presence of rotation, 
and in particular for systems rotating sufficiently fast. Rotating systems were 
tackled in the past by different models (see e.g. [11, 12] for a discussion of 
the stability of rotating gaseous cylinders, and the more recent works [13, 14] 
and references therein) , but a complete analysis of the Antonov problem with 
angular momentum is still lacking. Some progress was achieved recently for 
another class of systems where particles are assumed to obey a statistics that 
prohibits overlapping, the so-called Lynden-Bell statistics [15,16]. 

Indeed, we will show that the above picture (a)-(e) holds with minor modifica- 
tions for slowly rotating systems, that is when \L\ = L is sufficiently small. In 
particular, defining the reduced angular momentum \ = L/ \J RGN^ , we iden- 
tify an Antonov line e^X), with ec(0) = ~ —0.335 which plays in rotating 
systems the same role as the Antonov point plays in non-rotating ones. For A 
larger than a critical value Ac ^ 0.45, however, the collapse occurs not to a 
single dense cluster but to a double cluster (a "binary star" ) , signaling a break- 
ing of the (axial-)rotational symmetry of (1). It is remarkable that, strictly 
speaking, for A > Ac there is no Antonov limit: double cluster solutions exist 
at all energies, albeit being locally unstable (they are saddle points of entropy 
in the space of p's). The lower the energy, the more dense and point-hke the 
two clusters. We calculate the global phase diagram of the model in the (e. A) 
plane and analyze the thermodynamics by studying the phase transition that 
take place and the caloric curves. The equilibrium density profiles obtained 
in the different regions of the phase diagram show that, depending on the 
values of e and A, the formation of "rings" , "single stars" , "double stars" and 
"disk" -like structures is possible. 

From a broad statistical mechanics viewpoint, it is well known that self- 
gravitating systems, and more generally systems with long-range interaction 



3 



potentials^, represent a subtle and delicate issue [17]. For a recent critical 
discussion of the general statistical mechanics of non-trivial systems see [18]. 
The main problem is that the usual thermodynamic limit {N oo,V ^ 
oo,N/V finite) for systems with potentials such as (2) does not exist [19]. 
Their thermodynamics is evidently not extensive. Now extensive systems nec- 
essarily have positive specific heat, but this need not be true for non-extensive 
ones [6,20-24]. Hence the canonical and microcanonical description of non- 
extensive systems do not necessarily coincide (they do only in "regular" phases 
with positive specific heat, and in the dilute limit (A?" ^ oo, F — > oo, N/V^^^ 
finite) discussed at length in [25,26]). In the most general and physically inter- 
esting cases, self-gravitating systems should be studied in the microcanonical 
ensemble, a circumstance that obviously introduces several technical difficul- 
ties which are amplified by the introduction of rotation. 

The choice of Boltzmann statistics, which is made by most authors [4], to 
describe the behaviour of particles contributes to complicate the situation 
further. The basic assumption is that all particles can eventually occupy the 
same point/cell in space, i.e. overlapping between particles is admitted. With 
Boltzmann statistics, the system (1,2) has no ground state, hence the entropy 
is not bounded. This is the ultimate origin of Antonov's catastrophe, even 
(as we shall see) in rotating systems. Clearly, if one wants to describe a sys- 
tem where Newtonian gravity is the dominant macroscopic interaction then 
configurations with high particle density should be avoided, for at such short 
distances quantum effects and nuclear interactions become more important 
than gravity. One way around this problem is to introduce hard cores for par- 
ticles. Another substantially equivalent way is to use Lynden-Bell statistics [27] 
instead of Boltzmann's. 

We shall proceed as follows. In Sec. 2, we will expound our microcanonical 
mean-field theory, deriving the analogous of Eq. (4) for rotating systems. In 
Sec. 3 we will show results from the numerical solution of the new equation, 
and discuss in particular the phase diagram in the (e, A) plane and the caloric 
curves. A comparison between the Lynden-Bell scenario discussed in [15, 16] 
and the Boltzmann theory is presented in Sec. 4. Finally, we formulate our 
conclusions with some final remarks . 



2 Mean-field theory 

The theoretical analysis of Antonov's problem in a rotating framework can be 
carried out along the lines of the mean-field approximation described in [15,16], 

^ By which we mean potentials decaying with the interparticle distance r more 
slowly than r~^~^ with 6 > in D dimensions when r — oo. 
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where the equihbrium properties of a self-gravitating and rotating system were 
studied assuming Lynden-BeU statistics for particles ^ . The micro canonical 
entropy in our case is given by 



Sn = Sn{E,L) = log 



a 

m 



N 

/ 6{Hn - E) S{L - ^ X p,) Dr Dp 

■' i=l 



(5) 



where a is a constant that makes the argument of the logarithm dimensionless. 
The integral over momenta, namely 



r I N N 

K,L)= 5{K - -^.p!) ^{L - ^ ^« x Pr) Dp 

1=1 i=l 



(6) 



can be calculated via its Laplace transform (integrals resulting from insertion 
of the integral representation of the 5-function are at most of Gaussian type, 
hence trivial [15,16,28]) 



F^({r,}, s,L)= F^({r,}, K, L) e"^^ dK {^s > 0) 
Jo 

3N-3 

= (detI)-V2(^) ^ ^-^Lr^-^L 

where I = II({r'j}) denotes the inertia tensor, with elements {a,b — 1,2, 3) 

N 



(7) 
(8) 



Iab({ri}) = ^(r- - ri,ari,b) 
1=1 

and L^I~^L — J2a,b Lalab^b- Inversion of (8) gives 



(9) 



3N-3 
(27r)-T- 
rSW-Ss 







{K-\L^I-^lY-^ K>\L'^l-^L 

otherwise 



(10) 



hence one is left with 
Sn = log 



oA r [E - lL^r'L-<l>({r,})Y 



-Dr 



where A = (27r)^/r((3iV - 3)/2). If one is interested (as we are) in the 
large N limit, one can keep only the terms of dominating order in N in (11). 



^ A calculation similar to the one presented here (with Boltzmann statistics) was 
performed in [28], where however the resulting stationarity condition is studied in 
the context of the Thirring model [6] , a simpUfied model of a rotating self-gravitating 
gas. 
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This yields 



■>N 



log 



— j[E- ^L^r'L - ^{r,})]-Dr 



(12) 



and it remains to integrate over V 



N 



To this aim, we partition the spherical box V into K identical cells labeled 
by the positions of their centers {k = 1, . . . , K), and we let n{rk) denote 
the number of particles inside the cell at position r^, with I]^iT^(t'a;) = A^. 
In order to turn the integral into a sum over the cells, we adopt the following 
mean-field approximation for the potential $ and the components of the inertia 
tensor as functionals of the density profile p = p{r) — {K/V)n{r): 



Iab{{ri}) Iab[p] = j p{r) (r^Sab - ran) dr 



(13) 
(14) 



This choice automatically excludes correlations from the theory. The micro- 
canonical mean-field entropy can now be written as 



'at 



log 



aA 
M 



\ f 1 

-j[E--L^I-^L-m'^P\p\dp{r) 



(15) 



where P[p] is the probability to observe a density profile p. Now P[p] is just 
proportional to the number of ways in which our particles can be distributed 
into the K cells, and is thus determined by the statistical model by which one 
describes the behaviour of the particles. Assuming Boltzmann statistics, there 
is no spatial constraint for particles so that all of them can eventually occupy 
the same cell in V . This leads to 



P[p]« 



N\ 



n{ri)\- ■ ■n{rK) 



N\ exp 



NK 



j c(r) (log c(r) — l)dr 



(16) 



where in obtaining the last expression we used Stirling's approximation and 
defined the relative cell occupancy c(r) = n{r)/N — Vp{r)/{KN). Introduc- 
ing the dimensionless variable x — r/R and neglecting irrelevant constants, 
one finally arrives at 



S 



mf 
N 



log 



dc{x) 



(17) 



S[c] = -log 



1 



2 ^ ^ 



c{x){\ogc{x)-l)dx (18) 



where © = / c{x)dx — 

For N large though finite, S"^^ can be estimated by the steepest descent 
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method, which for S = S^^N yields 



5" = max 

c{x) 



log 



E--Ln-'L-^\c\ 
2 ^ ' 



1 

e 



c(a;) (log c(a;) — l)dx 



(19) 



Notice that the stationarity properties of -^jv do not depend on ©. This is most 
easily seen by setting f{x) — c{x)/Q, which leads to 



S — max 



3, 

2 log 



E-^L^r'L-m'^ - I f{x){\ogf{x)-l)dx + \oge 



(20) 

We remark that I = In terms of / and in units of GN'^/R and NR"^, 
respectively, $ and lab (see (13), (9)) become simply 



^■'^ 2 J \x-x'\ 



(21) 

(22) 



The stationarity condition 6S/6f = implies the following equation for the 
extrema of the entropy: 



f{x) = exp 



\x — x' 



dx' 



oj X xY + A* 



(23) 



where a; = I is the angular velocity and // is a Lagrange multiplier that 
ensures the constraint / f{x)dx — 1, and 



3/2 



[E - \Ln-^L - $[/]] 
is the inverse kinetic energy. Notice that (3 = 



(24) 



Eq. (23) is the analogous of (4) for a rotating system. Unfortunately, it cannot 
be solved analytically. However, some simple qualitative considerations are 
possible. In absence of rotation, the solution of (23) is known to depend only 
on X = \x\. The minimum expected effect of rotation is to create distortions 
that cannot be described by x alone. Hence it is necessary to introduce angular 
variables, and it is convenient to pass to spherical coordinates {x,6,(f)). If the 
system is slowly rotating and if, as we shall assume in the following, the angular 
momentum is parallel to the 3-axis, one may neglect the 0-dependence and 
write formally the solution of (23) as a sum of Legendre polynomials weighted 
by functions of x alone ^ : 



f(x) = J2M^)Pi{^o8e) 

1=0 



(25) 



^ When in the following we speak of axial rotational symmetry or for short of axial 
symmetry we refer always to the 3-axis. 
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This procedure is well-known [29]. and has been recently applied to the Antonov 
problem for slowly rotating systems by Chavanis [14]. The latter, however, do 
not display any serious qualitative difference from non rotating ones. The most 
interesting phenomena occur at high angular momenta, and in order to de- 
scribe them it is necessary to include the dependence on 0. We thus write the 
solution of (23) as 

oo I 

fix)=Jl E fl,mi^)yimiO,<P) (26) 
1=0 m=-l 

where Yi^ are real spherical harmonics and fi^jn are their weights (again func- 
tions of X alone), which will have to be determined. Clearly, this expansion 
differs from the one in Legcndre polynomials essentially for the terms with 
m 7^ 0. We shall see that these terms are the ones that provide the relevant 
new physics at high angular momenta. 

Using this expansion it is easy to show the following statements. The normal- 
ization / f{x)dx = 1 is equivalent to 

2^ j x'U,o{x)dx = 1 (27) 

so that the (0, 0) harmonics is ensures this condition. The coordinates of the 
center of mass are given by 

xf^ = j xj{x)dx = 2^1 x'f,,i{x)dx 

J X2f{x)dx = 2^11 x^h,-^{x)dx (28) 



x?^ = J xsf{x)dx ^'^sJlJ x^fiA^)dx 



so that harmonics with / = 1 are responsible for fixing the center of mass. It is 
convenient to set the latter in the center of the spherical box (a; = 0), which 
amounts to postulating /i „j = for m e {—1, 0, 1}. The diagonal components 
of the inertia matrix read ^ 



-^22 — 
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J x^fo,o{x)dx + J x'^f2,o{x)dx - \j J x'^f2,2{x)dx (29) 
■ j x'^fofl{x)dx + y ^ / x'^f2fl{x)dx + Y / ^'^h2i^)dx (30) 



-^33 — 



J x'^fo,Q{x)dx - J x^f2,Q{x)dx (31) 



^ We warn that indices like in I^b refer to spatial components, while indices like 
/; ^ refer to harmonics. 
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while for the off-diagonal elements one gets 




J x^h,-2{^)dx (32) 
g / x'hMdx (33) 
^ / x'h,.,{x)dx (34) 



In particular, one has 



-^11 — 122 



j x'^f2,2{x)dx (35) 



so that this parameter, which will play a foremost role in the following analysis, 
only depends on the (2, 2) harmonics. 

Using spherical harmonics as a basis set, for the Newtonian potential one 
has [30,31]: 

^ Ji = t E . ^^'^irnie', ct>') (36) 



\x-x'\ ^^tl;2/ + lmax{a;,x'}^+i 

Now it is easy to see that the so far unknown functions //^^ satisfy the following 
integral equation (which can be found by multiplying both sides of (23) by 
Yijn and integrating over angular variables): 

fl,ra{x) = j Yira{e, 0) e^T.V -I' M-)'^l'rn'{e,'t>)HP'-'-' ^ dOdct) (37) 

with 

47r /■ min|x,x'V „ . ,x9 , , , ^ 

^lAx) = 7^7— r / / n\+i iAx'){x'fdx' 38 

2/ + 1 J max{x,x'}'+^ 

Despite their unfriendly look, Eqs (37,38), which should be solved together 
in a self-consistent way, can be dealt with numerically. Technical details are 
similar to those described in [15, 16] in relation to the Lynden-Bell statistics 
case, and wc will skip them here. Let it suffice to say that for the computation 
the series (36) was truncated at the order /max = 12 and that all harmonics 
(even and odd) up to /max were included, except those with 1 = 1. Notice that 
this does not imply a truncation of the expansion (26), i.e. the fi^rn exist also 
for orders higher than /max- The effects of truncating the potential are well- 
known, and are discussed e.g. in [16]. In short, when truncated at higher and 
higher order the series (36) differs more and more from the real potential in 
a narrower and narrower region around x = x', in a way that is reminiscent 
of the Gibbs phenomenon. As we said, the angular momentum was taken to 
lie along the 3-direction. In this setup, solutions of (23) were calculated as 
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functions of e = and A = ^rgn^ ' ^^^^^^^ ^^^^ analysis are given and 
discussed in the next section. 



3 Results 



3.1 Relevant quantities 



The crucial quantity in the analysis that will follow is given by 

Hes(e, A) = det = K1K2 (39) 

.dedxS dlS 

where k,i^2 are the eigenvalues of the matrix, ordered assuming k,i > K2- The 
relevance of Hes(e, A) lies in the fact that it relates the thermostatistical equi- 
librium properties of systems in the microcanonical ensemble to the topology 
of the entropy surface in the space of conserved macroscopic quantities. In gen- 
eral, it can be shown [17,32] that pure, stable thermodynamic phases with pos- 
itive heat capacity are characterized by having Ki < (hence Hes(e, A) > 0). 
Here, 5 as a function of both e and A is concave, and statistical ensembles (mi- 
crocanonical and canonical) are equivalent. When at least ki > 0, one has a 
phase coexistence region with negative specific heat, and statistical ensembles 
are inequivalent. S{e, A) has a convex part that is ignored by the canonical en- 
semble, which considers the concave hull of S. The convex region corresponds 
to negative specific heat. For self-gravitating systems, (39) provides a power- 
ful characterization of some of the relevant quantities. In particular, we shall 
concentrate on the following aspects of the phase diagram. 

Antonov limit. For each fixed A, the Antonov limit ec(A) is the minimum e 
for which a solution of (23) exists. In general, ec(A) can be defined by the 
condition 

lim Ki = +00 (40) 

(Strictly speaking, Ki can diverge only in the mean-field approximation in the 
— > 00 limit at fixed e and A.) This is a straightforward generalization of the 
non- rotating case, where at the Antonov point Cc = ec(0) one has 

d3 

lim^ = -Foo (41) 

elec de 

When two (or more) conserved quantities have to be taken into consideration, 
the above condition becomes a condition on the maximum eigenvalue of the 
matrix of second derivatives of S with respect to the conserved quantities. 
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Upon varying A, the Antonov limit ec(A) describes a line in the (e, A) plane. 
We call this line the Antonov line. 

Phase coexistence. By virtue of the properties of (39), starting in the high- 
energy phase and decreasing e at fixed A one finds the beginning of a phase 
coexistence region when 

Hes(e, A) = (42) 

is fulfilled. We shall denote by eg (A) the corresponding value of e. In particular, 
at eo(A) the specific heat at constant (3uj = — |^ (a; = |a;|), i.e. 

''^'^ ~ Hes(e, A) ^A^ 

has a pole and changes sign (from positive to negative) upon decreasing e. 
Again, varying A the point eo(A) describes a line in the (e. A) plane (with 
eo(0) ~ —0.2), which separates a pure thermodynamic phase from a phase 
coexistence region. 

Bifurcation to axially asymmetric solutions. For sufficiently high A, solutions 
of (23) without axial rotational symmetry ("binary stars") bifurcate (con- 
tinuously) from axially-symmetric solutions [15,16]. This bifurcation can be 
detected with the aim of the order parameter (see (35)) 



D = |7n - 722 

introduced in [15]. It is easy to understand that solutions with D = Q are 
axially-symmetric, while solutions with D ^ {) break axial rotational symme- 
try. If A = the system is necessarily isotropic (/n = I22 = hs) and spherical 
symmetry cannot be broken. This reflects the fact that the solution of (23) 
depends only on x, so that the series expansion (26) is truncated after the first 
term, which involves only the Iqo harmonics. When A 7^ 0, anisotropics may 
occur (733 7^ 7ii,722) and one can have cither solutions that are rotationally- 
symmetric around the 3-axis (7ii = 722 or 7? = 0), which eventually depend 
on X and 9 to account for obvious rotation-induced distortions, or solutions 
that are asymmetric around the 3-axis (7ii 7^ 722 or 7? 7^ 0). The latter must 
depend on x, 9 and 0, and correspond to double clusters. However, it turns 
out that such solutions appear only when A is larger than a threshold value 
and e is sufficiently low (see also the discussion at the end of Sec. 4, and in 
particular Fig. 5). 

Stability of the double cluster solutions. Finally, we shall be interested in the 
local stability properties of double-cluster solutions (high-e, axially symmetric 
solutions are locally stable because in this limit the kinetic term of (1) domi- 
nates and the system behaves more and more as an ideal gas). This amounts 
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(43) 
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Fig. 1. Microcanonical phase diagram. GAS (axially symmetric) denotes the pure 
phase; phase coexistence regions with negative specific heat are denoted by PC-I 
(single-cluster plus "gas") and PC-II (double-cluster plus "gas"). Lines are as fol- 
lows. Solid-continuous: Antonov line ec(A), Eq. (40) . Light-dashed: phase boundary 
eo(A), Eq. (42). The continuous line from A to B marks the beginning of the bifurca- 
tion line for the onset of double cluster solutions and separates the PC-I and PC-II 
phases; this line merges with eo(A) right of A. Heavy-dashed: marginal stability line 
for the double cluster solution. Results have been obtained with Zmax = 12. 



to looking for the marginal stability line defined by the condition 



G[f*]=0 



(45) 



where 



G[/] = fix) - exp 



J \x — X 



dx' 



UJ X cc)^ -I- \x 



(46) 



and /* is an axially asymmetric solution of (23). The marginal stability line 
marks the border between the region of the phase diagram where axially asym- 
metric solution are stable from that where they are at most metastable (saddle 
points in the entropy surface). In principle, at the most general level the lo- 
cal stability analysis requires solving the eigenvalue problem for the kernel 
Kix, x') = TTT-iTTT-TT- We have not been able to carry out this task. However, 
it is possible to get numerical information about (45). Some technical details 
about this problem are given in [16], and we will not discuss them here. 
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3.2 Phase diagram 



The phase diagram resulting from numerical solution of (23) is displayed in 
Fig. 1. In Fig. 2 we show instead a small sample of equilibrium density profiles. 
The situation for sufficiently low angular A resembles that of non-rotating 
systems. At high e, a pure "gas" phase is found, where the equilibrium density 
profiles are axially-symmetric and particles are uniformly spread over V . As e 
is decreased, the Newtonian contribution becomes more and more important 
and at the eQ^X) line, the specific heat becomes negative and one enters a 
phase coexistence region, labeled PC-I, where a collapsed configuration (see 
Fig. 2a) competes with the gas state (see Fig. 2b). Upon further decreasing e 
no more solutions of (23) are found and the Antonov limit is reached. In this 
region, we indeed find that 

lim Hes(e, A) = -oo (47) 

4ec{A) 

which, comparing with (40), implies that /t2 < 0. Rotation has the only effects 
of shifting both the onset of negative heat capacity and the Antonov limit to 
higher e's (as also argued in [14]), and of deforming the density profile (as 
clear from Fig. 2b). 

When A becomes larger than the critical value Ac — 0.45, however, the picture 
changes drastically. Here axially-symmetric solutions become unstable against 
fluctuations that break such a rotational symmetry and axially-asymmetric so- 
lutions, i.e. "double stars", occur (see Fig. 2c). The picture is clear for high A. 
Upon decreasing e, one encounters the eo(A) line (which in that region coin- 
cides with the bifurcation line) and enters a phase coexistence region (labeled 
PC-II) with negative specific heat cp^i, Eq- (43)), where double-cluster and gas 
configurations compete. However, at odds with the low-A further de- 




Fig. 2. Example of stable equilibrium configurations, contour plots and (above) 
density profiles. Left to right: single-cluster (collapsed state, e = —0.33, A = 0.001); 
"disk" (e = -0.1, A = 0.4); double-cluster (e = 0, A = 0.6). 
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crease in e makes double-clusters unstable. Remarkably, the solution for lower 
e still exists, hence there is no Antonov limit for A > Ac strictly speaking. This 
issue will be further discussed in the next section. 

The bifurcation line and the eo(A) line do not coincide in a small range close 
to Ac (the AB line in Fig. 1 gives the bifurcation line in this range). Here 
the bifurcation line separates the two mixed phases, i.e. upon crossing it by 
increasing A at fixed e double clusters are obtained. Transitions from PC-I to 
PC-II can also occur upon decreasing e at fixed A. 

It is important to stress that these results were obtained by truncating the po- 
tential via (36) to order /max — 12. While calculations performed with smaller 
Imax in the lower half of the phase diagram (A < Ac) gave essentially the same 
results, for high angular momenta we observed some dependence on Imax, al- 
beit weak. In particular, it is possible that upon increasing Zmax the marginal 
stability line is further inclined towards the eo(A) line. As a consequence, the 
PC-II phase would become thinner and thinner, so that double stars would 
be stable (if at all) only in a very tiny portion of the parameter space. This 
quantitative issue, which is strictly related to the use of Boltzmann statis- 
tics, certainly deserves further investigation. We believe that no other main 
(qualitative) changes to the picture presented here should occur. 

Finally, moving on to thermodynamics, in Fig. 3 we show the caloric curves 
/3 vs e at fixed low and, respectively, high A. The peak in (b) corresponds to 
the bifurcation line. 



4 Compcirison with Lynden-Bell statistics 



It is interesting to compare the behaviour of the self-gravitating and rotat- 
ing gas of Boltzmann particles described above with that of a same system 
of Lynden Bell particles [15,16]. Using Lynden-Bell statistics [27] instead of 
Boltzmann's amounts to postulating a maximum capacity of no <^ iV parti- 
cles for each of the K cells in which V is subdivided in order to compute the 
integral (12), so that Eq. (16) would become 

pr 1 ^! TT ^o' 



n{ri)\ ■ ■■n{rK)\ J^^^ K - n{rk))\ 

From the physical viewpoint, Lynden-Bell systems are leptodermous. This is 

not possible for Boltzmann particles, as seen by the fact that the clusters 
in Fig. 2 are sharply peaked, with the density decreasing smoothly from the 
maximum (Lynden-Bell particles can form stars with a dense, finite core and a 
fast-decreasing density [15]). Correspondingly, the stationary point condition 
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Fig. 3. Caloric curves for (a) low and (b) high A. The inset in (b) shows the behaviour 
of the order parameter at A = 0.7. ebif in (b) denotes the bifurcation energy at 
A = 0.7. 
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for c = ©/ reads 

" P f <x') 1 



c{x) = 



1 + exp 



P f C[X' , , i ^, ,2 



This case has been studied thoroughly in [15, 16]. Notice that, at odds with 
what happens in Boltzmann statistics (which is the uq —>■ N followed by 
N ^ oo limit of the Lynden-Bell case), here the result depends on ©, which 
is now given by © = uoKB? ^ ^^"^ related to the total mass of the system. 

The micro canonical phase diagram for this case is discussed in [15, 16], where 
it is shown that the system has three pure phases ("gas", single cluster and 
double-cluster) separated by a large phase coexistence region with negative 
specific heat. For historical reasons, in this paper we shall discuss the canoni- 
cal phase diagram, which can be obtained by mapping the (e. A) variables onto 
their conjugate intensive variables (T = l/j3,uj). Evidently, in the canonical 
ensemble phase coexistence regions cannot be detected [17], as one must al- 
ways have positive specific heat. This means that the collapse temperature is 
determined by the value of (3 at which the derivative of (3 versus e changes sign, 
and the region PC-1, corresponding to lower energies, is not accessible. For the 
present comparison purposes, however, this limited information is sufficient. 
Canonical phase diagrams are shown in Fig. 4. 

In the Boltzmann case, one sees that for high temperatures the equilibrium 
state is an axially-symmetric, gas-like one. At a critical temperature Tduj) the 
system collapses into a single cluster (for a; < — 1-03) or into a double 
cluster (for uj > ujc) and lower temperatures are forbidden. (For a; = the 
known value T'c(O) ~ 0.3967 is recovered.) Lynden-Bell statistics turns out to 
shift the onset of both the collapse to a single cluster and the formation of 
the double cluster to lower temperatures. However, at lower temperatures the 
system can exist in a single-cluster state (low u) or in a double cluster state 
(high uj) . Notice how in both cases the phase coexistence regions with negative 
specific heat that are found in the microcanonical phase diagrams are lost in 
the canonical setup. 

Nevertheless, the behaviour of the Boltzmann system beyond the marginal 
stability line is quite intriguing. As we said, Eq. (23) admits double-cluster 
type of solutions for all e if A > Ac — 0.45, so there is no Antonov limit. 
However, such solutions are locally unstable: there exist perturbations of the 
density profile which lead to an increase of entropy. A profile of the entropy 
density S versus e at fixed A = 0.46 is shown in Fig. 5. At high energy the 
system is found in a stable, axially symmetric state: a "gas". Upon decreasing 
e, the specific heat CQ^^ (43) stays positive until point C is reached, where 
it turns negative (the curvature Ki of S'(e, A) or equivalently the quantity 
Hes(e, A) is zero). One enters the phase coexistence region where collapsed 
configurations are possible. A further (small) decrease in e will drive us to the 
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Fig. 4. Canonical phase diagram for (a) Boltzmann statistics and (b) Lyndcn-Bell 
statistics. The dashed line in (b) reproduces (a) and shows the difference between 
the two types of statistics. 
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Fig. 5. (Entropy profile as a function of e at fixed A = 0.46. Low energy branches 
correspond to double-cluster solution, the high-energy branch to the "gas" -like equi- 
libria. Notice that C is a maximum of S{e). 

bifurcation point D, where axial- rotational symmetry breaks down in favor of 
such structures as "binary stars". (The axially symmetric solution however 
persists, though locally unstable.) Decreasing e again one arrives at E, where 
double clusters become locally unstable, that is a saddle point of S[f]. Double 
cluster solutions (possibly more than one) however still exist left of E. In such 
states the two "stars" are extremely dense and almost point-like. 



5 Conclusion and outlook 

In summary, we have studied the equilibrium properties of a self-gravitating 
and rotating gas assuming that particles obey Boltzmann statistics in space. 
The physical properties of this system can be understood easily as deriving 
from the balance of (a) kinetic and gravitational energy at low angular momen- 
tum, and (b) gravitational and rotational energy at high angular momenta. 
In ), the effects of rotation are basically negligible except for trivial 

distortions of the equilibrium shapes. If the total energy is sufficiently high so 
that the kinetic term dominates the gravitational one, the equilibrium state of 
the system is axially-symmetric and gas-like. At low energies, instead, gravi- 
tation induces the system to collapse into a dense core in a low-density gas. In 
case (b), the high energy regime is similar but at low energy rotation is strong 
enough to hinder the collapse into a single cluster and drive the system to 
organize into a double cluster structure (a "binary star"). This simple, intu- 
itive picture is confirmed entirely by our mean-field theory, which shows that 
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in case (b) at least one main bifurcation (from axially symmetric to axially 
asymmetric) in the solutions of the basic Eq. (23) occurs, and possibly many 
others can be found upon varying the conserved quantities e and A. 

Our work leaves several directions for further research. The main question 
concerns probably the dynamics of these systems, which is an equally old [33] 
and widely-discussed theme in astrophysics, sec e.g. [34,35]. Among the many 
open issues and possible extensions for the equilibrium case (such as the effect 
of different masses, see [36,37]), we want to put forward the following few. 
First is the behaviour of the system left of the marginal stability line where 
double cluster solutions become unstable. (We mentioned in Sec. 3 that the 
location of this line turns out to depend on the maximum order of harmon- 
ics included in the calculation of the gravitational potential, /max- Also the 
fact that no Antonov hmit exists for A > Ac could be an effect of trunca- 
tion.) At low angular momenta, namely below Antonov limit, no statistical 
equilibrium state exists and dynamical methods should be used to grasp the 
physics of the system (see e.g. [5] and references therein). But for A > Ac it 
is still possible to obtain (unstable) equilibrium states as saddle point of the 
entropy surface. This means that certain perturbations applied to such "sin- 
gular" double-cluster structures (highly concentrated) will lead to an increase 
of entropy. It is not clear to us whether local entropy maxima (metastable 
states) exist in this region. Definitely, there are no global maxima for S is 
everywhere unbounded. This issue certainly deserves to be studied further. 
Another important question that we didn't explore concerns the relevance of 
mass as a further conserved quantity. If N is fixed, it should be treated on 
the same level as e and A, and this could lead to some qualitative changes in 
the picture presented here. As a final consideration, we stress the mean-field 
character of our solution. In spite of several claims being made of mean-field 
theory being exact for self-gravitating systems, we believe that it would be 
important to investigate the effects of particle-particle correlations, also at 
the merely quantitative level. 

Acknowledgments. Useful discussions with O. Fliegans are gratefully ac- 
knowledged. 
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